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Abstract 

We describe an exact derivation of the total nondissipative transverse force 
acting on a quantized vortex moving in a uniform background. The derivation 
is valid for neutral boson or fermion superfluids, provided the order parameter 
is a complex scalar quantity. The force is determined by the one-particle 
density matrix far away from the vortex core, and is found to be the Magnus 
force proportional to the superfluid density.^ 

Using the analogy to the classical fluid dynamics one finds that a mov- 
ing vortex in a superfluid feels a force proportional to its velocity vy as 
Fjvf = p s K x vy. H Here p s is the superfluid density and K the quantized 
circulation] 2 ] In earlier work we have shown the close connection between the 
Magnus force and a Berry phase 0, but our arguments might depend on 
a detailed model of the vortex structure. This does not give a satisfactory 
answer to the question of various possible corrections to this expression, as 
have been suggested in Ref. In the present paper we consider a single, 

isolated vortex forced to move through an infinite superfluid with uniform 
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velocity vy by some moving pinning potential. We are able to make an exact 
evaluation of the coefficient of vy in the force which has to be applied to the 
vortex in terms of the integral of the momentum density round a large loop 
that surrounds the vortex. The first step in the argument involves a pertur- 
bative calculation of the effect of the motion of the vortex in terms of the 
instantaneous eigenstates of the Hamiltonian. In the case of a system with no 
inhomogeneous substrate this expression can be written in terms of a Berry 
phase ||. This expression can be written in terms of the derivatives of one- 
and two-particle Dirac density matrices, or, alternatively, as the commutator 
of two components of the total momentum operator. This can be turned into 
an integral over a surface (or a loop for the two-dimensional case) surrounding 
the vortex and at a large distance from it. From this we can conclude that 
there are no corrections to p s JZ x vy. 

For simplicity we work with a two-dimensional system, but the argument 
can readily be applied to a three-dimensional system. We assume that the 
system is homogeneous and infinite. The isolated vortex in the superfluid 
can be pinned to a position ro by applying a potential J2i ^( r o — r i) which 
acts on all the particles in the superfluid. For a static problem this potential 
can be arbitrarily weak, but when the pinning potential is made to move 
through the superfluid it must have a strength sufficient that the vortex does 
not get detached from the pinning center, either by tunneling with the aid of 
the Magnus force, or by acquiring enough energy from the phonon system. 
The potential has to be strong enough to break the degeneracy of this broken 
symmetry vortex state and allow us to use a perturbative treatment of the 
velocity. The pinning potential may also be very strong, as it is in the Vinen 
experiment ||, where the circulation is pinned to a solid wire. Our results 
are quite independent of the details of the pinning potential, provided that 
it is strong enough to allow perturbative methods to be used, and does not 
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break cylindrical symmetry. 

In terms of the instantaneous eigenvalues E a {t) and eigenstates \^ a (t)), 
for which we choose phases such that (^ a \^a) = 0, and letting Vo denotes 
the partial derivative with respect to the position ro of the pinning potential, 
we have the expectation value of the force on the pinning potential first order 
in vortex velocity as 



iTlPn 

V H- ^-v v • Vo + h.c. 



Since Vo-ff is the commutator of the operator Vo with H, the commutator 
cancels the energy denominator, and eq. (1) can be rewritten as 



dxn 



3^ a 



dx 

where n is the unit vector normal to the plane. The Berry phase associated 
with a closed loop in ro can be written as the integral over the area enclosed 
by the loop of F/fivy. Equation (|2|) corresponds to the familiar form for the 
Berry phase ||. Since we can choose the wave functions in such a way that 
the dependence on ro is entirely through (r — ro), the partial derivatives with 
respect to ro can be replaced by a sum over partial derivatives with respect to 
the particle coordinates va. Upon thermal average, this expression can now 
be written in terms of the Dirac density matrices as 

Fjvy = -ihh ■ J d 2 r[V x Vp(r\r)] r=r > 

-ihh- J d 2 n J d 2 r 2 [2Vi xV' 2 r(ri,4;ri,r 2 )] r=r - , 

where p and T are the one- and two-particle Dirac density matrices for the 
system, and the sum over i and j denotes a sum over all the particles in the 
system. 



The integral over the two-particle density matrix T vanishes. A formal way 
of getting the result is to say that the first line of eq. (||) gives the transverse 
force in terms of the expectation value of the commutator of the x and y 
components of the total momentum of the particles. This commutator is a 
one-particle operator given by [P x ,P y ] = J J dxdy _ ^-|^J , where 

the ip<, ip are creation and annihilation operators for fermions or bosons. Since 
it is a one-particle operator, its expectation value is given by the one-particle 
Dirac density matrix. 

The integrand of the first term in eq. @ can be written as half the curl 
of (V — V)/?(r', r) evaluated at r = r'. Now divide the integrals over r and ri 
into a sum over finite areas labeled with an index a; the first of these areas 
will be centered on the position of the vortex, and the others are all well away 
from the vortex. Stokes' theorem can be used to write the result as 

This result is exact, and there are no contributions to the integrals from the 
neighborhood of the vortex core, since we have chosen the boundaries of the 
regions of integration to be well away from the core; any contributions from 
the core states or from the properties of the core must be reflected in the 
density matrices well away from the core. In particular, because localized 
states inside the vortex core, such as occur for an s-wave superconductor [|J, 
do not influence the one-particle density matrix far away from the core, their 
contributions to the transverse force must be zero. It should be pointed out 
that this property is not transparent in eq. (||), where both detailed forms of 
density matrices and explicit integrations are needed to recover it. 

For a neutral superfluid the integrand in eq. (^) is just the momentum 
density, equal, by definition, to p s u s + p n u n , where u s and u n are the local 
values of the superfluid and normal velocities. If the circulation of the normal 



fluid is zero, this term gives F = p s J£ x Vy, the Magnus force Fm Q^ll- 
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